The effective potential of the confinement order parameter in the Hamiltonian 

approach 
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The effective potential of the order parameter for confinement is calculated for SU(A^) Yang-Mills 
theory in the Hamiltonian approach. Compactifying one spatial dimension and using a background 
gauge fixing, this potential is obtained within a variational approach by minimizing the energy den- 
sity for given background field. In this formulation the inverse length of the compactified dimension 
represents the temperature. Using Gaussian trial wave functionals we establish an analytic rela- 
tion between the propagators in the background gauge at finite temperature and the corresponding 
zero-temperature propagators in Coulomb gauge. In the simplest truncation, neglecting the ghost 
and using the ultraviolet form of the gluon energy, we recover the Weiss potential. Neglecting the 
ghost and using for the gluon energy aj(p) the approximate Gribov formula aj(p) ~ p -|- M^/p one 
finds a critical temperature of ^fZMj-K. We explicitly show that the omission of the ghost drasti- 
cally increases the transition temperature. From the full non-perturbative potential (with the ghost 
included) we extract a critical temperature of the deconfinement phase transition of 269 MeV for 
the gauge group SU(2) and 283 MeV for SU(3). 

PACS numbers: ll.lO.Wx, 11.15.-q, 12.38.Aw, 12.38.Lg 



I. INTRODUCTION 

One of the major challenges of particle physics is 
the understanding of the deconfinement phase transition 
from the confined hadronic phase with chiral symme- 
try spontaneously broken to the deconfined quark gluon 
plasma with chiral symmetry restored. This transition 
is expected to be driven by the gluon dynamics, which 
suggests to investigate this transition first in pure Yang- 
Mills theory. In the absence of quarks the deconfine- 
ment phase transition is related to the center of the 
gauge group 1] : Center symmetry is realized in the low- 
temperature confining phase and spontaneously broken 
in the high-temperature deconfining phase. When quarks 
are included, center symmetry is explicitly broken and 
the deconfinement phase transition is expected to become 
a crossover. 

Understanding the deconfinement phase transition re- 
quires non-perturbative methods. In quenched QCD re- 
liable results have been obtained by means of the lattice 
simulations [3|- These methods become, however, ex- 
tremely expensive when dynamical quarks are included 
and fail at large baryon densities due to the notorious 
fermion sign problem at non-vanishing chemical poten- 
tial. Alternative non-perturbative methods, which are 
based on the continuum formulation of QCD are there- 
fore desirable. In recent years substantial progress has 
been made within continuum approaches to QCD [J-Q- 
Among these is a variational approach to the Hamilto- 
nian formulation of Yang-Mills theory in Coulomb gauge 
[6|-l8|. In this approach the energy density is minimized 
using Gaussian type ansatze for the Yang-Mills vacuum 
wave functional. Within the approach of Ref. [6] a de- 
cent description of the infrared sector of Yang-Mills the- 
ory was obtained |9l-[l4{ . Extension of this approach to 



full QCD has shown that the coupling of the quarks to 
the transversal gluons aniplifies the spontaneous break- 
ing of chiral symmetry [15|. Recently this approach 
was also extended to finite temperatures by considering 
the grand canonical ensemble, making a suitable quasi- 
particle ansatz for the density operator and minimizing 
the free energy [16l . [l7[ . In the present paper we present 
an alternative Hamiltonian approach to finite tempera- 
ture Yang-Mills theory, which does not require an ansatz 
for the density operator. The finite temperature is in- 
troduced here by compactifying one spatial dimension. 
From the physical point of view this approach there- 
fore addresses, strictly speaking, the Casimir problem. 
A summary of preliminary results obtained in this ap- 
proach was reported in [18] . 

The ground state properties of a quantum field the- 
ory are usually studied by means of the Euclidean space 
functional integral. This approach, which is also the ba- 
sis for lattice calculations, can be easily extended to finite 
temperatures by compactifying the Euclidean time axis, 
whence the momentum variable (energy) corresponding 
to the compactified Euclidean time becomes discrete with 
values given by multiples of lix times the temperature 
(Matsubara frequencies). Due to the 0(4) symmetry of 
Euclidean quantum field theory (at zero temperature) all 
four Euclidean dimensions are equivalent and instead of 
compactifying the time one can equally well compact- 
ify one of the spatial dimensions. This procedure can be 
also applied to the Hamiltonian formulation of a quantum 
field theory. Thus, instead of considering the canonical or 
grand canonical ensemble of a theory, its extension to fi- 
nite temperatures can be equally well obtained by consid- 
ering the theory for one compactified spatial dimension 
and identifying the length of the compactified dimension 
with the inverse temperature. In the present paper we 
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apply this procedure to study Yang-Mills theory at fi- 
nite temperature and, in particular, will investigate the 
deconfinement phase transition. To make contact with 
the previous Hamiltonian formulation of Yang-Mills the- 
ory in Coulomb gauge we use here a background gauge 
fixing, which reduces to Coulomb gauge for vanishing 
background field. We will consider a constant Abelian 
background field living in the Cartan subalgebra and be- 
ing directed along the compact dimension. We will then 
calculate the energy density as a function of the back- 
ground field and the length of the compactified dimen- 
sion. As will be shown in Sect. |ll] this quantity provides 
the finite-temperature effective potential of the Polyakov 
loop, which is the order parameter of confinement. From 
this potential the critical temperature of the deconfine- 
ment phase transition will be extracted. 

The organization of the paper is as follows: In the 
next section we show how the effective potential of the 
Polyakov loop can be calculated in the Hamiltonian ap- 
proach. In Sect. Ill Al we discuss how a background field 
chose in the Cartan subalgebra and directed along a com- 
pactified dimension can figure as an order parameter for 
confinement alternativ to the Polyakov loop. In Sect. Ill Bl 
we show how in the Hamiltonian approach a Lorentz in- 
variant quantum field theory can be extanded to finite 
temperature by compactifying a spatial dimension. The 
upshot of Sect.|lT]is that the effective potential of the con- 
finement order parameter can be obtained in the Hamil- 
tonian approach form the energy density calculated as 
a function of a constant background field living in the 
Cartan subalgebra and beeing directed along a compact- 
ified spatial dimension. The background field method 
is formulated in the Hamiltonian approach in Sect. IIIII 
In Sect. HVl we develop the Hamiltonian approach in the 
presence of a background field in a background gauge. 
The crucial point is here the resolution of Gauss' law. 
In the absence of the background field this approach re- 
duces to the Hamiltonian approach in Coulomb gauge 
[a|- We also discuss here the choice of the variational 
wave functional required to ensure that the expectation 
value of the gauge field agrees with the background field. 
The Dyson-Schwinger equations (DSE) for the various 
propagators in the presence of a background field are 
derived in Sect. |Vl In Sect. IVII we consider a constant 
background field living in the Cartan subalgebra of the 
gauge group and relate the solutions of the correspond- 
ing DSEs to the ones of the Hamiltonian approach in 
Coulomb gauge in the absence of the background field. 
For pedagogical reasons we develop the formalism for the 
gauge group SU(2) and present its extension to SU(Af) 
in Sect. IVIDI In Sect. IVIII we summarize the results ob- 
tained in the Hamiltonian approach in Coulomb gauge 
for the various propagators, which serve as input for the 
background field calculation of the effective potential of 
the Polyakov loop given in Sect. I Villi A short summary 
and our conclusions are given in Sect. IIXI Some lengthy 
mathematical derivations are deferred to appendices. 



II. THE EFFECTIVE POTENTIAL OF THE 

CONFINEMENT ORDER PARAMETER IN THE 

HAMILTONIAN APPROACH 

A. The Polyakov loop potential 

Consider Yang-Mills theory at finite temperatures and 
assume that the temperature is introduced by compact- 
ifying the euclidean time axis, identifying the length L 
of the compactified dimension with the inverse temper- 
ature. In the absence of fermions the finite-temperature 
Yang-Mills theory is invariant under gauge transforma- 
tions C/(xo,x), which are temporally periodic up to an 
element Zk of the center Z(iV) of the gauge group SU(A^) 

C/(i,x)=Zfe(x)C/(0,x), (1) 

where the center element is given by 

Zk = e-^''^"' = e'^'^l e Z(iV), fc = 0, 1, . . . , iV - 1 (2) 

with jlk — p-^Ta being a co-weight vector. (We use 
antihermitean generators Ta-) The order parameter of 
confinement is the expectation value {P[Ao]{x)) of the 
Polyakov loop J3} 



P[Ao](x) = — trT'cxp -/ dxoAo(xo,x) 



(3) 



This quantity can be related to the free energy -Fq(x) of 
an isolated static (infinitely heavy) quark at space posi- 
tion x by [l|Q 



(P[Ao](x))=c 



_ „-LF,(x) 



(4) 



and is entirely defined within pure Yang-Mills the- 
ory. Under gauge transformations satisfying Eq. ((T)) the 
Polyakov loop transforms as 



P[Ao](x)^ZfeP[Ao](x) 



(5) 



and so does (P[Ao]) -^ Zk{P[AQ]). It follows that the 
order parameter (P[Ao]) vanishes in the center symmetric 
confining phase, while (P[Ao]) 7^ in the deconfining 
phase with center symmetry broken. 

In the continuum theory the Polyakov loop is most 
easily calculated in Polyakov gauge, defined by OqAq — 



^ The phases of Yang-Mills theory can be alternatively distin- 
guised by the "disorder" parameter of confinement, the 't Hooft 
loop I20II . which is dual to the temporal Wilson loop but can be 
used both at zero and finite temperature. Its continuum repre- 
sentation was derived in Ref. [21]. Evaluation of the 't Hooft 
loop in the Hamiltonian approach in Coulomb gauge (at zero 
temperature) has given a perimenter law llll as expected for the 
confined phase. 



and Aq living in the Cartan subalgebrao In this gauge 
the path ordering in Eq. ([3]) becomes irrelevant 



B. Finite temperature from compactification of a 
spatial dimension 



P[^o](x)--tre-^^«W 



(6) 



and under a global center transformation ([SJ Aq trans- 
forms as 



Aq^ Aa + — /ife . 



(7) 



In Refs. [2J, [2^ it was shown that in this gauge the in- 
equality 



{P[Ao]) < P[{Ao)] 



(8) 



holds, and furthermore that also P[{Aq)] vanishes in the 
confining phase. This establishes P[{Ao)] (alternative to 
{P[Aq])) as an order parameter of confinement. Further- 
more, due to the unique relation ^ between the gauge 
field Ao in Polyakov gauge and P[Ao], instead of P[{Ao)] 
one can also use (Aq) itself as alternative order parame- 
ter of confinement. As a consequence the phase structure 
of finite temperature Yang-Mills theory can be extracted 
from the effective potential e[(Ao)] of (Aq). By gauge in- 
variance a non- vanishing (Aq) requires the presence of an 
external background field a. Choosing a to be constant 
and diagonal (thus being in the Polyakov gauge) and to 
satisfy a = (Aq) the background field a becomes an or- 
der parameter of confinement, whose value is determined 
by the minimum of the effective potential e[(Ao) = a]- 
Gauge invariance ensures that this potential is invariant 
under the Z{N) transformation ([7]) 



(^o> 



2n , 
1 



■fJ-k 



= e[{Ao) 



(9) 



This potential was first calculated in one loop perturba- 
tion theory in Ref. 26] . Corrections to this potential from 
the renormalization of the kinetic terms of the gauge field 
were calculated in Ref. [23|. It was found that e[a] has 
its minimum at a = 0, so that (P) ~ P[o] = 1, implying 
that the theory is in the deconfining phase with center 
symmetry broken. This is the expected behavior at high 
temperature, where the one-loop calculation is reliable. 
In the low temperature confined phase one expects that 
e[a] is minimal at "center symmetric configurations" a, 
for which (P) c^ P[a] = 0. Such a behavior was found e.g. 
in Ref. [23, l28[ , where the effective potential was deter- 
mined non-perturbatively from a renormalization group 
fiow equation. 



^ In Polyakov gauge topological defects like magnetic monopoles 
and center vortices occur 123 . 12311 . which presumably are the "con- 
finer" of the theory. 



Obviously, the order parameter (P) ~ P[{Ao)] or (Aq) 
is not directly accessible in Weyl gauge Aq =0, which is 
assumed in the canonical quantization. However, by 0(4) 
symmetry, all four Euclidean dimensions are equivalent 
and instead of compactifying the time, one can equally 
well introduce the temperature by compactifying one of 
the spatial dimension, say the xa-axis, and consider (A^) 
as an order parameter for confinement. This can be seen 
as follows: 

Consider Yang-Mills theory at finite temperature L^^, 
which is defined by the partition function 



Z(L)=Tre-^^(^) 



(10) 



Here H{A) is the usual Yang-Mills Hamiltonian defined 
by canonical quantization in Weyl gauge Aq = 0. The 
partition function (|10p can be equivalently represented 
by the Euclidean functional integral, see for example 
Ref. [H 



ZiL) 



a;*^'— pbc 



Y[VA^{x)e 



-S[A] 



where 



rL/2 p 

S[A]^ / dx^ / d3x£(A'^;a;^) 

J-L/2 J 



(11) 



(12) 



is the Euclidean action and the functional integration is 
performed over temporally periodic fields 



A^(L/2,x) = ^''(-L/2,x), 



(13) 



which is indicated in Eq. (|lll) by the subscript x'^ — pbc. 
This boundary condition is absolutely necessary at finite 
L but becomes irrelevant in the zero temperature {L — >■ 
oo) limit. 

We perform now the following change of variables 



a;0^z3 


A"^C3 


xi->zO 


A'^C° 


x^^z^ 


A^^C' 


x^-^z' 


A'^C' 



(14) 

Due to the 0(4) invariance of the Euclidean Lagrangian 
we have 

£(A'^,a;'')=/:(C^,z^) (15) 

and the partition function pip can be rewritten as 

Z{L)= I []2?C^(z)c-^[^"l, (16) 

where the action is now given by 

r rL/2 

SIC] = dz" dz^ dz^ dz^C{C^',z'') (17) 

J J-L/2 



and the functional integration runs over fields satisfying 
periodic boundary condition in the z'^— direction 



C^ (0", z\ z2, L/2) = C^ (z°, z\ z2, -L/2) 



(18) 



We can now interprete z° as time and z = (z^, z^, z^) as 
space coordinates and perform a usual canonical quan- 
tization in "Weyl gauge" C° = 0, interpreting C — 
{C^ ,0"^ ,C^^ as spatial coordinates of the gauge field, 
which are defined, however, not on R'^ but instead on 
IR-^ X S^ . We obtain then the usual Yang-Mills Hamil- 
tonian in which, however, the integration over z^ is re- 
stricted to the interval [—■§,■§■]■ Let us denote this 
Hamiltonian by H{C,L). Obviously H{C,L — > 00) = 
H{C). Reversing the steps which lead from pUj) to pT|) 
and taking into account the irrelevance of the temporal 
boundary conditions in the functional integral for an infi- 
nite time-interval we obtain from Eq. (|16p the alternative 
representation of the partition function 

Z(L) = Tr e- ^ d.«H(c,L) ^ ^-^ ^r e'^^^^'^) . (19) 

T— ^00 

Due to the infinite 2;''-(time-)interval T — > 00 only the 
lowest eigenvalue of H{C, L) contributes to the partition 
function Z{L). The calculation of Z{L) is thus reduced to 
solving the Schrodinger equation H{C, i)V'(C) — Eijj{C) 
for the vacuum state on the space manifold R^ x S^{L), 
where S^{L) is a circle with circumference L. 

The upshot of the above consideration is that finite- 
temperature gauge theory can be described in the Hamil- 
tonian approach by compactifying a spatial dimension 
and solving the corresponding Schrodinger equation for 
the vacuum sector. As the above derivation shows this 
equivalence holds for any 0(4) invariant quantum field 
theory. 

The above consideration for the partition function can 
be extended to the finite-temperature effective potential 
e[(Ao)]. One finds that e[{Ao) = a] can be calculated in 
the Hamiltonian approach from e[(A3) = a] with the z^- 
axis compactificd. Furthermore, as shown in Ref. |29| . in 
the Hamilton approach the effective potential eKA^) — a] 
is given by the energy density in the state that minimizes 
(H) for given (A3). In this paper we calculate the effec- 
tive potential e [(A3) = a] in the Hamiltonian approach 
using the representation (J19p of the partition function. 
In the spirit of the discussion given in Sect. Ill Al we will 
refer to £[(^3)] as the effective potential of the Polyakov 
loop. 

A few comments are in order: In the present paper we 
use the Hamiltonian approach merely as a mathematical 
tool to calculate the effective potential of the background 
field at finite temperatures. As the previous consider- 
ations show, strictly speaking, time in the present ap- 
proach is not the true time but a spatial dimension and 
the compactified spatial dimensions represents the time. 
Alternatively, interpreting the time of the present Hamil- 
tonian approach as the true time the present approach 
actually describes Yang-Mills theory in a space with one 



compactified spatial dimension and is thus reminiscent of 
the Casimir problem. 



III. THE BACKGROUND FIELD METHOD IN 
THE HAMILTONIAN APPROACH 

Let ^[A] = {A\ip) be the gauge invariant wave func- 
tional of the Yang-Mills vacuum. The generating func- 
tional of the Hamiltonian approach is then defined by 



sW'W = (V'|e^J^|V'). 



(20) 



Fixing the gauge by the standard Faddeev-Popov 
method we have explicitly 

e^y] = fvABetM[A]S{r[A])e-^^^^+^^^. (21) 



Here we have rewritten the wave functional as 



V'[A]=exp(--5[A] 



and assumed that it is real. Furthermore 

r[A]-o 

is the gauge fixing constraint and 

,at_sr[A^] 



M' 



se'' 



(22) 



(23) 



(24) 



is the corresponding Faddeev-Popov kernel with 

A'^ = UHd + A)U = U'' DU , C/ = exp(-e) (25) 

being the gauge transformation. Formally Eq. (PT|) looks 
like the standard Yang-Mills functional integral of the 
Lagrangian formulation in Euclidean space. However, 
in Eq. (PT|) S[A] is not the usual classical action and 
the functional integral runs over the spatial compo- 
nents of the gauge field only. In principle, S[A] is not 
known. However, in recent years there has been substan- 
tial progress in determining this wave functional approxi- 
mately from a variational principle. Furthermore, in 2+1 
dimension more detailed information on the Yang-Mills 
vacuum wave functional has been obtained from lattice 
studies [201 ^^d strong coupling considerations [3l|, [32] ■ 
In the present paper we will use for S[A] the results of the 
variational calculations in Coulomb gauge carried out in 
Refs. [a,B[l3- For the moment let us assume that ip[A] 
and thus S[A] is known. The similarity of Eq. (PT|) with 
the standard Lagrangian Yang-Mills functional integral 
allows us to apply the standard background field method 
|33| to simplify the evaluation of the effective action 



T[A] 



-W[j] 



JA 



A 



J=3iA) 



SW[j] 



(26) 



For this purpose we split the gauge field A into a back- 
ground field a and a fluctuating part A 



A^a + A. 



(27) 



The gauge transformation ([^5)) of the total gauge field A 
can be distributed among the background field a and the 
fluctuation A in different ways. The two common and 
convenient possibilities are: 

1. quantum gauge transformations, which leave the 
background field invariant 

SqQ = 

SgA = [D, 6Q] = [d + A, (56] , (28) 

where 

d = d + a (29) 

is the covariant derivative with the background field 

2. background gauge transformations 



Sbd — [d, Q] 
5bA = [A, 6] 



(30) 



Following the standard background field method [33J we 
introduce the generating functional 



^W[a,j] 



VA Det M 



[A]s[nA]) 



'S[a+A]+fjA. 



(31) 

which explicitly depends on the background field a since 
the source couples only to the fiuctuation A and since 
the gauge fixing functional /""[A] may depend on the 
background field. The action S[A] = S[a + A] is, of 
course, invariant under both transformations (|28p . (1301) . 
It is, however, convenient to eliminate the gauge degrees 
of freedom by fixing the invariance with respect to the 
quantum gauge transformation (pS)) . Then the Faddeev- 
Popov operator is given by 



M'^'iA] 






(32) 



As is well known, the generating functional Vl^[a, j] is in- 
variant under the background gauge transformations pop 
if the gauge fixing functional /"[^J transforms covari- 
antly under these transformations. A convenient gauge, 
which satisfies this requirement, is the background gauge 

r[^] = [d,A] = [d,A-a]=0, (33) 

for which the Faddeev-Popov operator p2p reads 



M[A] 



-Dd. 



(34) 



where the hat " " " denotes the adjoint representation of 
the gauge group and D is the covariant derivative 



D 



ah ra6 



S^'d + A 



ab 



A 



ab 



fab^c 



T: 



ab 



I 



act 



(35) 



with /"^"^ beeing the structure constants. This gauge is 
the 3-dimensional spatial analogue of the Landau-deWitt 
gauge. It is not difficult to see that in the background 
gauge ([55]) both the gauge constraint and the correspond- 
ing Faddeev-Popov determinant are invariant under the 
background transformation ([50]) . Consequently the gen- 
erating functional l/l^[a,j'] (l3T]) and hence also the effec- 
tive action 



5W[a,j] 



T[a,A]^ [W[a,j]- I Ai] .A 



Sj 



are invariant under the background transformation Sb 
pop. Furthermore, from the translation invariance of 
the integration measure VA follows that the true effec- 
tive action r[A], see Eq. ([25)1 . is related to r[a,^] by 



T[a,A] ^T[a + A] 
and the classical fields are related by 
A = a + A. 



(37) 
(38) 



Obviously the invariance of Ffa, A] under the background 
transformation ([50]) ensures the gauge invariance of r[A]. 
In this paper we will determine the effective potential 
r[a] = r[a,y^ = 0] for a constant background field a. 



IV. HAMILTONIAN APPROACH IN THE 
PRESENCE OF AN BACKGROUND FIELD 

Below we develop the Hamilton approach to Yang- 
Mills theory in the presence of an external background 
field. This approach is an extension of the Hamiltonian 
formulation in Coulomb gauge [a, |3J| and reduces to the 
latter when the background field is switched off. Al- 
though in the presence of the background field the Hamil- 
tonian approach could, in principle, be also formulated in 
Coulomb gauge the use of the more general background 
gauge ([55)) turns out to be advantageous. 



A. Generalities 

To simplify the bookkeeping we will use the compact 
notation A^^(xi) = ^(1) for colored Lorentz vectors like 
the gauge potential and an analogous notation for colored 
Lorentz scalars like the ghost C"-^{xi) = C(l). We also 
define in coordinate space 



5(1,2)-5'^1'^M(X1-X2) 



(39) 



A repeated label means summation over the discrete color 
(and Lorentz) indices along with integration over the d 
spatial coordinates 

d N^-l 

A{l)B{l) ^ d^xJ2T. ^^W^fcW- (40) 

-' k=l a=l 



Furthermore, indices will be suppressed when they can 
be easily restored from the context. 

The usual canonical quantization of gauge theory as- 
sumes Weyl gauge Aq = Q and results in a Hamiltonian 



H[A] 



g^n' 



1 



B[Af 



(41) 



where 11(1) = —i6/6A{l) is the momentum operator 
and B[A] is the non-Abelian magnetic field. Here we 
have absorbed the coupling constant g in the gauge 
field, gA — > A. This Hamiltonian is invariant under 
time-independent gauge transformations. In Weyl gauge 
Gauss' law escapes the (Heisenberg) equation of motion 
and has to be imposed as a constraint on the wave func- 
tional 



t>U^[A] = p,^t^[A] . 



(42) 



Here poxt is the charge density of external matter fields. 
Since DII is the generator of (time-independent) gauge 
transformations, in the absence of matter fields the wave 
functional has to be gauge invariant V'i^'^] = V"!^]- 

We are interested here in the Yang-Mills vacuum in the 
presence of an external vector field a. Let {ipa) denote the 
vacuum wave functional in the presence of the external 
field a and 

(...), :=(V'a|...|^a) (43) 

the corresponding vacuum expectation value. We as- 
sume here that the states jV'a) are properly normalized, 
(V'a|'0a) = 1- In the absence of external fields the gauge 
invariance of the vacuum wave functional |V'o) guaran- 
tees that the expectation value of all gauge dependent 
quantities vanishes. This refers, in particular, to the ex- 
pectation value of the gauge field 



(^01^1^-0)= 0. 



(44) 



The background field assigns a prescribed expectation 
value to the gauge field. If ipo[A] = {A\ipo) denotes 
the gauge invariant vacuum wave functional in the ab- 
sence of the background field (which necessarily satisfies 
Eq. pi)) ). then it is trivial to see by a change of variables 
that the wave functional 

MA]-=MA-a-] (45) 

has the property that 

{A)^^a. (46) 

Although the wave functional (P5|) satisfies (|46p , in gen- 
eral it will not minimize the energy under the constraint 
(|46p . The wave functional jV'a) which minimizes the en- 
ergy under the constraint (j46p 

(i/),^min, {A)^^a (47) 

can be obtained from the constrained variational princi- 
ple 

(i7-j(l)A(l))^^extr, (48) 

where j is a Lagrange multiplier chosen such that Eq. (j46p 
is fulfilled. 



B. The background gauge fixed Hamiltonian 

Instead of working with gauge invariant wave func- 
tionals, it is usually more convenient to fix the gauge, 
where those gauges are preferable which allow for an ex- 
plicit resolution of Gauss' law. The price one pays for 
the gauge fixing is that the Hamiltonian acquires a more 
complicated, usually non-local form. A convenient gauge 
in this respect is Coulomb gauge [3J|. However, in the 
present paper, where we are interested in the energy den- 
sity of the vacuum in the presence of an external back- 
ground field a. Coulomb gauge is not the optimal choice. 
In this case it is more convenient to choose the back- 
ground gauge ([55)1 which gauge fixes the fluctuating field 
A — A — a with respect to the background field a. 

In view of the gauge (I33p it is appropriate to intro- 
duce generalized longitudinal and transversal projection 
operators 



lij{x) ^ d.i (dd 
where 



j , t^j 



1 ^o,b cab c 

(49) 



df = S^'d, 



= r 



(50) 



is the covariant derivative (1341) with the total gauge field 
replaced by the background field. For a constant back- 
ground field, these projectors have the same properties 
as the ordinary longitudinal and transversal projectors 
as far as their spatial indices are concerned. However, 
they are non-trivial matrices in color space. Using these 
projectors we split the gauge field and the momentum 
operator into "longitudinal" and "transversal" part^ 



H = ni 



H^, 



^i — ^ijAj 



H 



liji-i-j ■ 



(51) 



The longitudinal part of the gauge field A'' will be later 
eliminated by the background gauge fixing, Eq. ([551) . 
while the longitudinal part of the momentum operator 
H'I is eliminated by resolving Gauss' law (|42)) . which can 
be explicitly done in the gauge p3p as we will show now. 
Inserting Eq. ([^T|) into Gauss' law P^ and solving for 
the longitudinal part of the momentum operator we find 



H"V 



:(-Dd) 



P-Ip: 



where 



P = Pcxt +Pdyn[^] 



(52) 



(53) 



^ We keep here the terms "longitudinal" and "transversal" al- 
though the components A" and A^ have this property only for 
a vanishing background field a = 0. 



is the total color charge density, which contains besides 
the external charge pext also the dynamical charge of the 
gauge bosons in the background gauge p3)) 

Pdyn[A] = -Dn^. (54) 

Rewriting the covariant derivative as 

D^d+(A-a] (55) 

and using dll = 0, the dynamical charge becomes 

Pdyn[^] = -(A-a)n^. (56) 

It depends only on the fluctuation A^A—a of the gauge 
field A around the background field a. Since the I, t (|49)) 
are orthogonal projectors we have 



Choice of the wave functional 



n^ = n 



n^ 



(57) 



Using this relation and Eq. ([5^ we find for the gauge 
fixed Yang-Mills Hamiltonian 

H[A] =1 (g^J-^[A]n^J[A]n-^ + ^B[A^]B[A^' 
■^ \ 9" 

+ Hc[A], (58) 

where 

J[A] = V)eiM[A] = Det(-Dd) (59) 

is the Faddeev-Popov determinant (cf. Eq. ([Ml)) and 

Iic[A] = ^^J-\A]p{l)J[A]F[A]{ia)pm (60) 

arises from the elimination of the longitudinal component 
of the momentum operator n" by means of Gauss' law, 
see Eq. (|5^ . and describes the interaction between the 
color charges. Here 



F[A] 



-Dd) i-dd] f-Ddl 



(61) 



is the analogue of the so-called Coulomb kernel [6j in 
the background gauge (j33p . For a vanishing background 
field a — Q the covariant background derivative d — d + a 
becomes the V-operator and the kernel (pT|) reduces to 
the ordinary Coulomb kernel [(y\. 

In the gauge-fixed theory the matrix elements of an 
observable 0[n, A] are defined by 



\o\ 



VASif[A])J[A]r[A] 0[n,A] ^[A] , (62) 



where f[A] — is the gauge fixing constraint ([55)1 and 
J[A] (Eq. (l59l) ) the corresponding Faddeev-Popov deter- 
minant (j59p . 



In Ref. Q a variational determination of the Yang- 
Mills vacuum wave functional was carried out (in the ab- 
sence of an external background field) in Coulomb gauge 
using the Gaussian type wave functionao 



MA] 
■0[A] := A/'exp 



J-'/'[AWA] 
1 
25 



-—A{IM1,2)A{2) 



(63) 



Here A/" is a normalization constant and uj is the varia- 
tional kernel. The ansatz (j63p has the advantage that it 
removes the Faddeev-Popov determinant J[A] from the 
integration measure (j62[) . Here we extend the variational 
calculations to the presence of an external background 
field a. To fuUfill the constraint (A)^ = a we use the 
trial wave functional 



MA] = j-'^^[A]m~a] 



(64) 



with iplA] defined in Eq. (p5)) . In the absence of the back- 
ground field Tpa [A] reduces to the vacuum wave functional 
"00 [A] (j63p . This trial wave functional satisfies already 
the constraint (|46)) . so it remains to minimize the energy 
(H)^ with respect to the kernel a;(l, 2). Due to the pres- 
ence of the colored background field the resulting kernel 
a;(l, 2) will be a non-trivial matrix in color space. 

With the trial wave functional (j64l) we find for the 
expectation value of any observable Op, A] from Eq. ([5^ 
after the shifi[3 {A — a) -^ A oi integration variables 



(O), ^ {MO[Tl,A]\^a) 



(65) 



VA6{f[A + a])iP*[A]0[n,A + a] iP[A] 
=:(0[n,A + a])o, 
where we have introduced the abbreviation 



d[u,A] = ji/2[^]0[n,A]j-i/2[^] 



(66) 



Equation (j65p is the vacuum expectation value of the 
observable (!)[H, A-f a] in the gauge fixed theory with the 
gauge condition, Eq. (|55| 



f[A + a]=dA = 0. 



(67) 



Note that in (O)^, ((S5|) . the gauge constraint ([57)) is im- 
plemented, which eliminates the longitudinal gauge field 
A' I defined by Eq. ([5T|l . Equation ([55]) applies in partic- 
ular to the gauge fixed Hamiltonian ([55)) . Note also that 



^ To be more precise in Ref. If| the ansatz I I63I I was used with J [A] 
II59II replaced by the Faddeev— Popov determinant in Coulomb 
gauge Det(— D9). Furthermore in the present case u](l,2) can 
have a non-trivial color structure. 

^ Note that the momentum operator El remains unchanged under 
the shift of coordinates [A — a) —^ A. 



after the shift of variables {A — a) 
of the gauge field (|5S|) becomes 



A the color charge which for the trial wave functional (1551) is given by 



Pdyn[A 



An^ 



(68) 



which is formally the same expression as obtained in 
Coulomb gauge [y, |3J| except that transversality is now 
defined by the gauge condition (p7| and thus depends on 
the background field a. 

In passing we notice that in the gauge fixed Hamilto- 
nian i/[n,A] ([55)) the gauge field enters only in form 
of the covariant derivative D = d + A. (Recall that 
the non-Abelian magnetic field can be written as Bi = 
^eijk[Dj,Dk]-) Therefore, after the shift of variables 
A ^ A + a the Hamiltonian H[Il,A + a] depends on 
the background field a only in the combination with the 
V-operator d + a = d. This is, in particular, true for the 
Faddeev-Popov determinant p4p 



J[A + a] = Det 



d + A d 



(69) 



and thus also for the transformed Hamiltonian iJ[n, A + 
a] defined by Eq. ([S5)) . 

We are eventually interested in the effective potential 
of the order parameter of confinement. For this purpose 
it is sufficient to consider constant background fields a, 
which we will assume from now on. Note, since a- a = 0, 
for constant background fields a we have d • a = and 
the gauge condition ([55)) reduces to the condition ([57)) 
d • A = 0. For the time being the color structure of the 
background field a — a''(x)T(, is arbitrary. 



V. YANG-MILLS DYNAMICS IN THE 
PRESENCE OF A BACKGROUND FIELD 

The presence of the background field influences the 
form of the propagators. In addition, the propagators 
depend on the gauge chosen. In the following we present 
the equation of motion for the ghost and gluon propaga- 
tors following from the constrained variational principle 
([47| in the background gauge (|67p with the trial wave 
functional ([M)) . Thereby we will make no assumption 
on the form of background field (except that it is con- 
stant). The resulting equations of motion are direct gen- 
eralizations of the equations obtained in the variational 
approach in Coulomb gauge 6] and reduce to the latter 
for a vanishing background field. 



A. The gluon propagator 



Since by Eq. ([65)) all expectation values (. . .)a in the 
state \ipa), satisfying the constraint ([^ . can be reduced 
to the vacuum expectation value (. . .)o it suffices to con- 
sider the gluon propagator 



I?(l,2) = (A(l)A(2))o 



(70) 



V{1,2)^!L.^-\1,2) 



(71) 



Independent of the form of the background field from the 
definition ()70p follows the symmetry relation 2? (1,2) — 
2?(2,l),i.e. 



I?,"f(x,y)=P?,"(y,x) 



(72) 



Since the gauge field is "transverse" (w.r.t. the covariant 
background derivative d ([50)) ) 

A^(x)=t^f(x)Af(x) (73) 

the gluon propagator is also "transverse" satisfying 



,bb' 



•\a'h' 



i^^,(x)t-(y)I?,^l',(x,y)=P-(x,y) 



(74) 



Note, the ordering is here important since i(x) is a differ- 
ential operator. By Wick's theorem expectation values in 
the state -00 [^] ([63)) can be entirely expressed in terms of 
the gluon propagator (1701) . see Sect. lVE)) 



B. The ghost DSE 

In the presence of the external background field the 
ghost propagator is defined by 



G^{M-M])a: 



(75) 



where A^[A] is the Faddeev-Popov kernel ([5^ and {■ ■ ■)a 
is defined in Eq. ([55)) . from which follows that the ghost 
propagator can be expressed as vacuum expectation value 



G^(M-^\A + a] 



d + A d 



(76) 



In the way described in Ref. !6?] one derives for this ghost 
propagator the following Dyson-Schwinger equation 



G-i=Goi-I], 



(77) 



where Gq is the bare ghost propagator, which is defined 
by Eq. ([75)) with ^ = 



Go = (^-ddj . (78) 

Furthermore, the ghost self-energy is given by 

I](l, 2) = ro(l, 3; 4)G(3, 3')r(3', 2; 4')I?(4', 4) , (79) 
where I?(l, 2) is the gluon propagator ([70)) and 

5M[A + a][l,2) 



ro(l,2;3) 



M(3) 



(80) 



is the bare ghost-gluon vertex. The full ghost-gluon ver- 
tex F is defined here by 



M-^\A + a]TQM-^\A + a]) = GFG . (81) 

/ 



Following Ref. we use the rainbow ladder approxi- 
mation replacing the full ghost-gluon vertex by the bare 
one. The justification for this is the following: As shown 
in Landau gauge 35] the ghost-gluon vertex is not renor- 
malized. This holds also in Coulomb gauge }10J. Further- 
more, in Coulomb gauge the dressing of the ghost-gluon 
vertex, at least to one- loop order, is small [1J|. We as- 
sume that this remains valid for the present gauge (1571) . 
which reduces to Coulomb gauge in the absence of the 
background field. 

With the explicit form of the Faddeev-Popov kernel 
M [A] dM]) we find for the bare ghost-gluon vertex (IMl) 

rg,fe(xi,X2;x3) 

= (fb4"(xi)5(xi - Xg)) di(xi)(5(xi - X2) , (82) 

with T^^ = /"'"=. Here we have suppressed the adjoint 
color indices ai , 02 , of the ghost legs and set 

r(l,2;3)^(r,3fc3(xi,X2;x3))'^^'^^ . (83) 

Replacing the full ghost-gluon vertex F by the bare one 
(15^ we find for the ghost self-energy ([7^ 

I](xi,X2) = [tli,{-^Mh^:i)DU (X3,xi)" 
xf, (4(xi)G(xi,X3))fb(di(x3)(5(x3,X2)) , (84) 
which is still a matrix in adjoint color space. 

C. The curvature 

The kinetic part of the transformed Yang-Mills Hamil- 
tonian H cf. (|55)) contains functional derivatives of the 
Faddeev-Popov determinant. In the vacuum expectation 
value of H (see subsection IV Dp these derivatives enter 
in form of the ghost loop 



X(l,2)- 



1 / 5'^\nJ[A + a] 



(85) 



2 \ 5A{l)6A{2) /o ' 

which in the present context is referred to as curvature 
Q. The curvature is defined here as in Ref. Q, however, 
with the argument of the Faddeev-Popov determinant 
shifted by the background field. With the definition of 
the ghost-gluon vertex ([5T|) it can be expressed as (cf. 
Ref. 6] for more details) 



x(l,2) = iTr(GF(l)Gro(2)) 



(86) 



Using again the bare ghost-gluon vertex approximation 
one finds with 



x^?(yi,y2)- \tli{yi)tw{y2) 



X tr 



Ta' 4'(yi)G(yi,y2) rfc-dp(y2)G(y2,yi) 



(87) 



where the trace is over the adjoint color space. 



D. The energy density and the gap equation 

Since the vacuum wave functional ip[A\ (|63p is Gaus- 
sian it is straightforward to calculate the energy {H) ^ 
(cf. (|65|) ) by using Wick's theorem and expressing 
{H [n, A -1- a])(j as a functional of the gluon propagator 
T) (|70D . For the Abelian part of the magnetic energy one 
finds 

1 



(^B)„ = ^[(-ddj(l,2)P(2,l') 
while the non- Abelian part gives 



1 



/ ttNA\ -"- rabc f: 

\^B )a= ^J f 



ah' c 



d'^a 



-,6b' 



T^ii {XjXjV^^^^XjX) 



■^b'c 



^b'cf 



+ V^^ [x, x)Vt^ {x, x) + Vf^ [x, x)Vt^{x, x)\ . (89) 

For the expectation value of the kinetic part Hk of the 
transformed Hamiltonian (1661) 



HK[A] = '-{n{i)n{i) 



1 / 6 



2 \5A{1) 5A{1) 
l5hiJ[A]5\nJ[A 
'4 5A{1) 5A{1) 



In J[A] 



(90) 



one finds, using the definition of the curvature ([55]) (see 
Refs. [y, [131 for more details) 



{Hk), 



P-1(1,1) 



-x(l,l) + x(l,2)P(2,3)x(3,l) 



(91) 



Variation of the energy {H) ^ with respect to the gluon 
propagator V (|70)) yields the gap equation. As shown 
in Ref. [ITI the Coulomb term Hq has little influence 
on the gluon sector. Furthermore, the non- Abelian part 
of the magnetic energy gives rise to a tadpole, which in 
the absence of the background field contributes a (UV- 
diverging) constant to the gap equation, which can be 
absorbed into a renormalization constant. Ignoring the 
Coulomb term and the non- Abelian part of the magnetic 
energy the gap equation arising from 



{HK + Hi)^^m\n 



(92) 



becomes 
9' 



■V-\2, 3)I?-^(3, 1) = i-tddt){2, l)+g\{2, 3)x(3, 1) . 

(93) 
Here we have also ignored the derivatives (5x/^2?, which 
give rise to two-loop terms since x itself contains already 
one ghost loop. Using the gap equation ([M]) to express 
(— dd) in the magnetic energy in terms of V and x the 
vacuum energy can be cast into the compact form 



{Hk + H-^)a = 9^ 



Vi(l,l)-x(l,l) 



(94) 
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In Sect. IVIIII we will extract from this expression the 
effective potential of the order parameter of confinement. 



VI. THE PROPAGATORS IN THE PRESENCE 
OF THE BACKGROUND FIELD 

The above considerations are valid for any gauge group 
and for arbitrary constant background fields. For ped- 
agogical reason we will confine ourselves below to the 
gauge group SU(2). The extension to SU(A^) is straight- 
forward and will be presented in Subsect. IVIDI 



A. Choice of the background field 

As shown in Sect.HJto find the effective potential of the 
confinement order parameter, it is sufficient to consider 
constant background fields a — const. By isotropy of 
space, without loss of generality we can direct the back- 
ground field along the 3- axis 



a ~ ae^ . 



(95) 



Furthermore, as we have seen in Sect. |TT] the background 
field has to live in the Cartan subalgebra to figure as 
order parameter of confinement. Therefore we have to 
choose for SU(2) 



a = a Tb — aT^ , a — const . 



(96) 



In the Hamiltonian approach presented above the back- 
ground field occurs in the adjoint representation. The 
Cartan generator in the adjoint representation T^ = e'^^^ 
is diagonalized in the basis of the spin s = 1 eigenstates 

f3\a)=-ia\a), a = 0,±l, (97) 

where, in the usual bracket notation, the 

{a\a) = el (98) 

are the cartesian components of the spherical unit vectors 
(in color space) 



1 



eo-=i 



V2 



1 , e„=_i 




Here and in the following Greek letters a,T, p, . . . denote 
spherical color components (1,0,-1) while Latin letters 
a,b,c, . . . denote the Cartesian color components (1, 2, 3). 
The vectors e^- satisfy the symmetry relation 



e; = (-)^e_. 
and form an orthonormal 3-bein 



(100) 



(101) 



By Eq. ^ we have 

a\a) — aT^\a) — ~ia<j\a) 



(102) 



Any matrix in the adjoint representation can be ex- 
pressed in the spherical basis as 



M"^ ^^{a\(j)Ar^T\h) . 



(103) 



Due to Eq. (|102p . the covariant derivative ([50)) dk be- 
comes diagonal in the spherical basis 

d^ = 5''^dl, dl=du+ 5uz{-\aa) . (104) 

The same is true for the longitudinal and transversal pro- 
jectors (gH) 



tki=^''''tkh tki=Ski-lki, (105) 



which makes the spherical basis favorable. Defining the 
momentum representation of these quantities by 



d{x)e'P'' = d(p)e'P^ 



(106) 



we have in the spherical basis 



dkip) = i (P - <^ae3),^ = i (Pfe - cra^g) , (107) 
dUP)dnP) 



tli{p)=S, 



kl 



d-(p)d-(p) 



(108) 



where 



- d'^(p)d"(p) = (p - aae^f = pi + (pg - oaf (109) 

and pj^ is the projection of p on the 1-2-plane. Further- 
more these quantities satisfy the symmetry relations 

d^(-p) = -dr(p)- (110) 



B. Propagators 

For a constant background field the homogeneity of 
space is preserved and the Green's functions can depend 
only on coordinate differences and thus can be Fourier 
transformed as in the absence of the background field. 
We define the Fourier transform for the gluon propagator 
by 



2?(x,y) 



dpe 



ip(x-y) 



2?(p), dp ; 



dp 

(27r)° 



(111) 



and analogously for the other two-point functions like 
ghost propagator and curvature. The background field 
(1961) . however, singles out a direction in 3-space and thus 
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spoils SO (3) invariance. Therefore the Fourier trans- 
formed two-point functions, hke 2?(p), will not just de- 
pend on the modulus |p| but also on the direction of p. 
Since the background field a = aT^ is diagonal in the 
spherical basis and since d is the only color dependent 
quantity we expect that the various two-point functions 
like gluon and ghost propagators are also diagonal in this 
basis. We show in appendix \K\ that the equations of mo- 
tion can indeed be consistently solved for color diagonal 
propagators of the form 



C(p) = (5""G"(p). 



(112) 
(113) 



For the reduced propagators G"^ (p) , V^ (p) the equations 
of motion simplify drastically. From the appendix |2] we 
find for the ghost DSE (|MT| 



G'^ip)-' = G^oiPr'-^''ip): 



(114) 



where 



G^(p)-i = -d^(p)d-(p) = (p - e.aaf (115) 

is the bare (inverse) ghost propagator ((78)) and 

E-(p) = 
-Y, fdqdr{p)t^,{q)dUp)'D^{cDG'^+^{p + q) (116) 

is the ghost self-energy (I84p . Here and in the following 
sums of spherical indices, like /i-l-cr, are defined modulo 3. 
As shown in the appendix El we find for the gap equation 
(jA33|) after the transversal projectors are contracted 



^ (V^ip))-' - -d-(p)d-(p) + / {x"{p)f , (117) 



where (cf. Eq. ([MT 



C. Relating the propagators in the presence of the 

background field to the vacuum propagators in 

Coulomb gauge 

The solutions of the constraint variational problem 
(1171) are defined by the coupled equations (J114p . (J117L 
in which the expressions (|116l) S'^(p) and (|118l) x'^(p) 
enter. It is not difficult to see that solutions exist, which 
obey the symmetry relations 

P''(-p)=2?-''(p), G^(-p) = G-''(p). (120) 

Indeed assuming that these relations hold, one finds from 
the explicit expressions for E (|116p and x (|118l) by using 

dnni) 

E^(-p) = I]-^(p), X'^(-P)=X"''(P), (121) 

which in turn entails Eqs. (|120p . We will now show that 
these solutions can be related to the vacuum solution in 
Coulomb gauge obtained in Ref. Q- 

Fourier transforming Eq. d7T|) . going to the spherical 
basis (I103|) and assuming (|112D we have 



V^ip) = 



2uj-{p) 



(122) 



where uj" (p) is the energy of a gluon with spherical color 
component a. Using also Eqs. (|115l) and Eq. (jA28p from 
appendix El the gap equation (|117l) becomes 



iLo^p)f = {p-aae,f + {x^p))\ 



(123) 



which is now a scalar equation. It differs from the vac- 
uum gap equation in Coulomb gauge obtained in Ref. Q 
by the explicit presence of background field a in the first 
term of the r.h.sO and by the index a labeling the eigen- 
values of — iTa. Obviously, the gap equation (I123P allows 
for solution of the form 



a;^(p)=w(pj 



(124) 



where 



X'^(P) = 



2{d-l) 
^y'd<7<(q)5fe-,(p)<(q)G'^(q)G-+^(p + q) (118) 



is the scalar curvature. Finally from Eq. (|A35|) we have 
for the energy (|94p 



{Hk + Hg ) a — 



P^ 



aae^ 



(125) 



and tt'(p) is the solution in the absence of the background 
field, i.e. for a = 0, provided also the curvature has the 
property: 



x"(p)-x(pj. 



(126) 



We will now explicitly show that both the gap equation 
(fng)) and the ghost DSE (fTTi)) are indeed satisfied by 
propagators of the form 



((i-l)F^y'dg (^^2?-(q)-i-x'^(q)). (119) 2?-(p)=P(pJ, G'^(p) = G(pJ 



(127) 



Let us emphasize that G'^(p), E'^(p), V^{p) and x'^(p) 
are all scalars in IR'^ but have a non-trivial color structure 
due to their dependence on the spherical color index a. 



As, compared to Ref. [0j], in the present paper we have ignored 
the tadpole and the Coulomb term. 
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where T>{p) and G(p) are the corresponding propagators 
for a = 0, i.e. the vacuum solutions of the Hamihonian 
approach in Coulomb gauge [y|. 

First notice that by the definitions (|107l) , (llOSp we have 

d'^(p) = rf(pj, i-(p)=t(pj. (128) 

Assuming (|127p and noticing that by Eq. (|125p 



(p + ^)a+^, = p^ + q^ 

from Eq. (fTT6)) follows indeed with ([728]) 
I]-(p) ^ E(pJ 



(129) 



(130) 



and thus the ghost DSE (I114p allows indeed for solutions 
of the form (jl27p . In the same way one shows that, as- 
suming ([T?f)) . from Eq. (ITT^ follows 



X"(P)=X(P.): 



(131) 



which implies that also the gap equation (J123I) allows for 
solutions of the form (|127p . This completes the proof of 
existence of solutions of the form (|127l) . Since the prop- 
agators in the absence of the background field depend 
only on p^ the solutions (I127P also satisfy the symmetry 
relation (|120p . We have thus shown that the propagators 
in the background gauge (|67p in the state of minimal en- 
ergy under the constraint (|i5)) are related by Eq. 
to the ordinary vacuum propagators in Coulomb gaug^ 
The latter have been determined previously in a varia- 
tional approach in the continuum theory (y, |9| as well 
as measured on the lattice [361 . In the section FVlIIi we 
will use these propagators as input to calculate the effec- 
tive potential of the confinement order parameter. For 



this purpose we briefly review in Scct. lVIII the results ob- 
tained for w(p) and x(p) ii^ the Hamiltonian approach in 
Coulomb gauge at zero temperature. 



D. Extension to SU(Ar) 

The previous considerations can be straightforwardly 
extended to an arbitrary gauge group SU(A^). As ex- 
plained in Sect.HJthe background field a has to be chosen 
in the Cartan subalgebra 



E«'^^' 



fc=i 



where 



Hk, k = l,. 



(132) 



(133) 



are the generators of the Cartan subalgebra with r being 
the rank of the group. Since the Hk commute with each 
other they can be simultaneously diagonalized. Their 
eigenvalues are the so-called weights — i/ife, which are the 
entries of the weight vectors 



M = (MI:---:/^!-) • 



(134) 



In the present approach we need the background field in 
the adjoint representation a. The weights in the adjoint 
representation Hk are the roots ak 



Hk\(j) = -i(Jk\cr) , 
which are the entries of the root vectors 

(T — {ai,. . . jCTr) • 



(135) 



(136) 



For the gauge group SU(2), which has rank r = 1, the 
single generator of the Cartan subalgebra is given by 
Hi = Ta and the eigenbasis of T3 is explicitly given in 
Eq. ([M]). As in the SU(2) case, in the Cartan basis (fHS)) 
the background field (|132p is obviously diagonal 



a\a) = — ia<T|(T) 
where a = (oi, . . . , a^) and 

r 

acr ^ y^flfcafc 

fe=i 



(137) 



(138) 



is the scalar product in the Cartan subalgebra. Further- 
more for SU(A^) the non-vanishing roots come in pairs 
±cr. (Some of the roots may vanish as in the SU(2) case. 
In the mathematical literature the term "root" is usually 
reserved to the non-vanishing roots.) Besides SU(2), we 
will be mainly interested in the gauge group SU(3), which 
has rank 2. The two generators of the Cartan subalgebra 
are usually chosen as Hi = T3 and H2 = T^. For this 
group the positive weights read 



fi 



"•;^ 



1 1 



1 



1 



X2V3J' V2' 2^/3 
and the non-vanishing positivqfl roots are given by 



(139) 



(1,0), 



1 1 
2' 2 



V3 



U 



(140) 



where the first root occurs also in the SU(2) case. Note 
that for SU(3) the scalar product of a root vector with a 
weight vector is either or 2 



cTfie{0, -}. 



(141) 



''' In principle, there could, of course, be other solutions which do 
not have the form II127I I. So far we have no indication that such 
solutions do exist. 



^ A weight or a root is called "positive" if its first component is 
positive. 
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All considerations of the previous section remain valid 
for SU(3) with the only modification that the shifted mo- 
mentum (|125|) has to be replaced by 



((Ta)e3 



(142) 



where era is defined in Eq. (|138p . With this modifica- 
tion all equations of the previous section remain valid 
when the index a — 0, ±1 is replaced in SU(3) by all 
(2-dimensional) root vectors a. Since the roots (|140l) are 
normalized crcr — 1 the co- weight vectors jlk (Eq. dH)) 
are given by 2/x and accordingly the non-trivial center 
elements read 



1000 



100 ■ 




0.001 



1000 10000 



-iiifiii 



(143) 



where Hi = T3 , H2 = Tg are the generators of the Car- 
tan algebra in the fundamental representation, normal- 
ized to tr(TaT;,) = —^dab- Center symmetry of the effec- 
tive potential ^ requires that it remains invariant under 
the replacement of the background field 



An 
a^a+— /x, 



(144) 



where /x represents one of the roots (|139p . 



FIG. 1. The gluon energy a;(p) and modulus of the curvature 
x{p) resulting from the full numerical solution of the varia- 
tional approach in Coulomb gauge as described in Ref. fl7| . 
The Figure contains also the Gribov formula (|151[) fitted to 
the numerical data for ijj{p), as well as the IR form xm (|152|) 
of the curvature. In addition we also show the parameteriza- 
tion H154|) Xpara(p) to the numerical results for x{p)y together 
with the solution a;para(p) of the gap equation (|145ll assuming 
the optimal parameterization Xpara(p) p54[) as input. 



with p = |p|. From the ghost DSE pl7)l one finds the 
sum rule 



VII. THE COULOMB GAUGE PROPAGATORS 



Since p^^g = p (see Eq. ([125)) ) by Eq. ([127)) the zero- 
temperature propagators in Coulomb gauge (for vanish- 
ing background field) are given by the background gauge 
propagators with a = 0. Putting cr = in Eq. ()123p 
the zero-temperature gap equation in Coulomb gauge be- 
comes in the presently used approximation 



^2(p) 



X'(P) 



(145) 



Since x(p) — X'^^°(p) (jUSp is defined in terms of the 
ghost propagator G'^(p) this equation has to be solved 
simultaneously with the ghost DSE ()114p . Introducing 
the ghost form factor d{p) by 



G^=^ip) = 



5P^ 



(146) 



the ghost DSE becomes 



1 



d-'{p) = -+Id{p), 



(147) 



hip) = N, J aq[l - ipqf] 



djp - q) 1 

(p - q)2 2L.(q) • 



Equations ()145p and ()147p can be solved analytically in 
the IR using the power law ansatze [a, [13, flTf 



uj{p) 



A 

pa 



dip) 



B 



(148) 



a = 2;9 - rf + 2 



(149) 



For d ^ 3 the gap equation ()145p allows then for two 
solutions with ghost IR exponents 



/3=1, 



/? = 0.795. 



(150) 



The first one, for which a = 1, is compatible with the 
gluon propagator measured on the lattice. Indeed, the 
lattice results [33 for the gluon energy a;(p) can be well 
fitted by the Gribov formula [32 



(151) 




^oip) = \ P^ + ^r 



with a mass M ~ 880 MeV for SU(2) [sg and about 
the same value for SU(3) [S^. Using for a;(p) the Gribov 
formula ()15ip we find from the gap equation ()145p for 
the curvature 



xip) = -■ Xmip) , 

P 



(152) 



which represents the correct infrared behavior of xip) 
obtained in Refs. [a, [ifl [ill- What is missing in the 
expression (|152l) is the UV-part 



Xuv(p) 



yjp^ + A p^oo^ p 

In P^+^ In P^ 



c = const , 



(153) 

which arises in a self-consistent solution of the coupled 
ghost DSE pi7)) and gap equation ([TI5)) [i,[l3- In the 
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gap equation (I145p this term is UV-subleading compared 
to the UV- leading term ui{p) ^ p. 

Alternative to the use of the lattice result we solve 
the ghost DSE (|147l) and the gap equation (|145p numer- 
ically in the way described in Ref. [1^], where also the 
renormalization of these equations is discussed in detail. 
The resulting numerical solutions for w(p) — w°'^'^(p) and 
x{p) = x'^^"(p) ^-re shown in Fig. [T]for the (3 — 1 solu- 
tion. The physical scale is determined as in Ref. [17] by 
fitting the numerical results for uj{p) to the Gribov for- 
mula (fTSTj) . using the lattice result of Af = 880 MeV as 
input. Fig. [T] shows also the fit of the the Gribov formula 
(|151l) to the numerical results for u}{p) and xir(p) (USD)- 
In the double logarithmic plot the Gribov formula (J15ip 
is indistinguishable from the numerical results. 

For practical reason in the evaluation of the effective 
potential of the background field we parameterize the 
numerical results for xip) by 



Xpara(p) = u{p)xm{p) + w(p)XUv(p) 



(154) 



where xir{p) and xuv(p) are defined by Eqs. p52p and 
(|153D . and u{p) and v{p) are smooth cut-off functions 
to restrict xmip) and xuv(p) to the IR and UV-regime, 
respectively. An optimal fit to the numerical results for 
xip) is obtained with the choice 



u{p) 



v{p) 



A 



p2 + A 



(155) 

From the best fit of Eq. (|154l) to the numerical data 
one extracts the following values for the parameters (for 
SU(2)) 



77 = 4.83M, A = 89.62M2 
m = 3.13, 71 = 0.89. 



(156) 



Fig. [T] shows also the fitted curve (|154l) together with the 
numerical data for x(p)- With w(p) and x{p) at hand we 
are now in a position to evaluate the effective potential 
of the Polyakov loop, see Eq. (|162p below. 



VIII. THE EFFECTIVE POTENTIAL OF THE 
POLYAKOV LOOP 

As discussed in Sect.|ll]a constant background field a 
living in the Cartan algebra can serve as order parameter 
for confinement when it is directed along a compactified 
dimension. Thereby the inverse of the length L of the 
compactified dimension figures as temperature. In the 
Hamiltonian approach the effective potential of a spatial 
background field a is given by the energy density [2!|] 



e(a, L) 



{H)a 



Vid-1) 



(157) 



where (H)^ is defined by the constrained variational 
problem (j47p and given in the present approach by 



Eq. (|119p . Furthermore, V is the spatial volume and for 
convenience we have divided the energy by the number 
of transversal spatial dimensions. 

Assuming d = 3 and compactifying the 3-axis to a cir- 
cle with circumference L the above formulae all remain 
valid except for the following modifications: The integra- 
tion over coordinate space reads 

y"dV(x)=y"dvy" dx3fi^^,x3), (i58) 

where xj_ denotes the projection of x onto the 1-2-plane. 
Accordingly the integration over momentum space be- 
comes 

Jd'pf{p) = Jd^P±j J2 /(P±'P")' (159) 



n— — 00 



where 



Pn 



2im 



(160) 



are the Matsubara frequencies and the shifted momen- 
tum variable p^, (|142p reads 



Va=V± + iPn - c^a) es : 



(161) 



With Eq. (jl59p we find from Eq. (|119p for the energy 
density (IT57| . using (fT22)) . (fT24)) and (fT26)) 

e(a,L) = ^ /'d3p(c.-(p)-x'^(p)) 

where pcr — |Pg.| is defined by Eq. (|16ip . By shifting the 
summation index n one easily verifies that by Eq. (J14ip 
the potential (jl62p indeed has the periodicity property 



e (a + An/Lfi, L) = e(a, L) 



(163) 



required for the potential of the Polyakov loop by center 
symmetry. 

For later use we calculate first the contributions to the 
effective potential arising from a gluon energy uj{pa) that 
obeys a power law 



.„(p,A)-Afi-"(Vp2 + A)" 



(164) 



For later convenience we have introduced here an ad- 
dition parameter A > 0, which of course is irrelevant 
for P(j — > 00. Furthermore, M is a constant of dimen- 
sion mass, which was introduced for dimensional reasons. 
For the /3 = 1 solution of the variational calculation (see 
Sect. IVIip the IR and UV behavior of a;(fc) is given by 
uia=-i{p, A = 0) and u)a=i{p, A = 0), respectively. 

The contribution of uja{paT^) to the energy density 



Saia,L,X) := E / <J^P-L-^^w„(p^,A) 



(165) 
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FIG. 2. The effective potential a{a,L) p62p multiplied by L* at different temperatures L~^ for the gauge group SU(2). The 
potential is shown as function of the dimensionless variable x H178|) . where acr — a^ = a for SU(2), and for various forms 
of the gluon energy a;(p) and neglecting the ghost loop, x(p) = 0- i^) UV-(Weiss) potential (|170|) ijj{p) = p, (b) massive 
dispersion relation H173P uj{p) = i/m^ + p'-^, (c) IR-potential (|176p uj{p) = h'P /p, (d) the sum of the UV and IR potential, 
io{p) = p + M^ Ip, which is obtained by adding the two potentials shown in (a) and (c). 



is calculated in appendix [B] One finds after the value of 
the potential at vanishing background field is subtracted, 



see Eq. (jBT2|) . 

ea{a, i. A) := ea(a, L, A) - 6^(0, i. A) 



(47r)3/2r I 



'HH 



nL 



sm 



""fi)A-_,_,(„Xv^), 



(166) 



where Kr{z) is the modified Bessel function (JBIII) . For 
A = the above expression simplifies to (see Eq. (jBlOp ) 



ea{a,L, A = 0) 



8M 



/fl — a 



(47r)3/2 r(-f) 



3+a 



oo _. 



ncrSiL 



^' ' ^- ) , (167) 



cr>On=l 

where the summation is now over the positive roots only. 



Before we give a full numerical calculation of the effec- 
tive potential e(a,i) p62p let us discuss its qualitative 
features to reveal how the deconfinement phase transition 
shows up in this quantity. For this purpose we calculate 
below the IR and UV contributions to the effective po- 
tential. Note also that the momentum integrals (I162p are 
UV divergent, so that the UV contributions require an 
analytic treatment to remove the divergences. 



A. Asymptotic contributions 

For a first qualitative discussion of the Polyakov loop 
potential let us ignore the curvature x(p) in Eq. (|162p . 
x{p) is UV-subleading but has the same IR behavior as 
a;(p), see Sect. I VIII Therefore, we expect a substantial 
change of the deconfinement phase transition tempera- 
ture by ignoring x(p)- Nevertheless, in order to see how 
the deconfinement phase transition manifests itself in the 
effective potential it is enlightening to consider the sim- 
plified potential arising from (jl62p when the curvature 
x(p) is ignored. We first calculate the effective potential 
arising from the UV behavior of w(p). 
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In the ultraviolet the gluon energy is given by uj{p) = 
IpI, which follows from Eq. (jl64p with a = 1 and A = 0. 
For these values we obtain from Eqs. (|165p . (|167p 



e'i^y{a,L) := ea=i{a,L,X = 0) -eQ=i(0, L,A = 0) 
8 1 Y^ ■^ sin^ (nacri/2) 



7r2L4 



(168) 



(t>0ti=1 



The summation over the Matsubara frequencies can be 
carried out for < a(TL/2TT < 1 using 



n=l 



which yields 



TT 

90 



4 71^ 



12 



Tra; 
T2" 



a; 
48' 



< x < 27r 
(169) 



"UV 



(«'^)=oTl 



3L4 ^^ 

cr>0 



acri 

2tt 



acrL 

~2^ 



-1 



(170) 



This is precisely the Weiss potential, which is shown in 
Fig. 2(a) for the gauge group SU(2). It was originally 



obtained in Ref. 26] in a 1-loop calculation in Landau 
gauge. Thus, the present quasi-particle type approxima- 
tion (jl65p to the gluons (cf. also Eq. ([S5]) ') with the gluon 
energy uj{p) replaced by its ultraviolet limit, i.e. by the 
photon energy uJa^i(j),X = 0) = |p|, is equivalent to an 
ordinary 1-loop background calculation. Since up to the 
periodicity ^ the Weiss potential is minimal only for 
a = we find for the Polyakov loop 



{P[A,])^P[{A,)^a^O] = l 



(171) 



indicating the deconfined phase. 

For the sake of illustration we also consider a massive 
gluon dispersion relation 



ujip)^^AP+p\ 



(172) 



which is obtained in the variational approach to Yang- 
Mills theory in Coulomb gauge above the critical tem- 
perature of the deconfinement phase transition [13] ■ This 
dispersion relation follows from Eq. (I164p for a = 1 and 
A = M^. With these parameter values we find from 
Eq. (jl66p for the effective potential 

ei(a, L, M^) := ei(a, L, M^) - ei(0, L, M^) 

/I #2 °" / 9 \ 2 

= ^EEb7 sin2(n^aL/2)i^_2(n£M). 



27r2 



nL 



(173) 



This potential is show in Fig. |2(b)| for SU(2). For M = 
this potential reduc es, of course, to the Weiss potential 
shown in Fig. 2(a) Due to the presence of the energy 



scale M in Eg. (jl72]) this potential is now temperature 
dependent. ( The W eiss potential L'^e'^y is independent 



of L, see Fig. 2(a) ) For any value of M this potential 





(b) 



FIG. 3. (a) The Weiss potential (flTO)) and (b) the infrared 
potential (fT76)) at L~^ = 290 MeV for the gauge group SU(3) 
as a function of the dimensionless fields variables (|178|l x, y 
and multiplied by L*. In both plots the local minima are 
marked by contours in the sy-plane. 



is minimal at a = and gives rise to the same Polyakov 
loop as the Weiss potential. 

The quasi-particle approximation to the energy den- 
sity is, however, a priori non-perturbative. The gluon 
energy tj(p) determined in the variational approach, in 
particular, captures the infrared confinement properties 
of the gluons. In the infrared the gluon energy (I15ip is 
given by Eq. (I164p with a = — 1 and A = 0. For these 
values we find from Eq. (|167p 

e5^(a, L) :— ea=-i{a, L, A = 0) — eQ,=_i(0, L, A = 0) 

JVf2 1 °° 1 

= "VZ^ EE;^ «-'(--- V2). (174) 

cr>0 n=l 

The summation over n can be carried out for < 
aaL /2ti < 1 using 



ECOS(nx) TT TTX X ^ „ /,„„x 

^^=6-T + T' 0^-^2., (175) 



n=l 



n^ 
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which yields 



AJ 
Er^{a,L)^2—Y: 

cr>0 



acrL 

2n 



BLcrL 

2tt 



(176) 



This potential is shown in Fig. 2(c) for the gauge group 
SU(2). It drastically differs from the Weiss potential 



(jl68p obtained from the UV behavior oiuj{p) and shown 
in Fig. 2(a) First it has the opposite sign and second a 



non-trivial L-dependence. Its minimum occurs at a = ^ 
corresponding to a center symmetric ground state. Ac- 
cordingly it yields a vanishing expectation value of the 
Polyakovloop (IT7T|) . 



{P[A3 



P 



{A,)^a^- 



= 0. 



(177) 



In Figs. [3] (a) and (b) the asymptotic UV and IR po- 
tentials (|170p and (|176p . respectively, are shown for the 
gauge group SU(3) as function of the dimensionless vari- 
ables 



X : = 



27r 



y 



2tt 



(178) 



By the periodicity of the effective potential these vari- 
ables can be restricted to the intervals 



0<x<l, 0<y< 



V3- 



(179) 



As mentioned before, the Weiss potential L'^e'^^{a,L) is 
independent of the temperature L~^, while the IR poten- 
tial is temperature dependent and is shown in Fig. |3(b)| 
for L^^ = 290 MeV. The minima of the Weiss potential 
are degenerate and occur at field configurations amin, for 
which the Polyakov line yields a center element z 



-.-La 



^ z € Z{N) , 



(180) 



and thus P[aniiiJ = 2 7^ as expected for the deconfined 
phase. For SU(3) the minima of the Weiss potential occur 
at (see Fig. 3(a)) 



(XmiTi, ymin) = (0, 0) , ( 0, ^= ) , I 1, ^= 



V3 



V3 



for which the Polyakov line is 



= 1, e'~l, e"'~ 



(181) 



(182) 



The minima of the confining IR potential occur at center 
symmetric configurations, for which the Polyakov loop 
vanishes. For SU(2) this configuration is given by Xmin = 
1/2 (see Fig. 2(c)[) while for SU(3) these configurations 
read (see Fig. 3(b) I: 
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FIG. 4. The ag = slice of the potential p87|) for the gauge 
group SU(3) as a function of x (|178p for different L~^ in the 
range from 484 MeV to 880 MeV (from bottom to top). The 
critical temperature is 597 MeV. 



and the corresponding Polyakov lines are given by 
e-La.„.„ ^diagfe'T", e-'^", 1 

diag f e"'T^, 1, e'^ 
diagfe'^, 1, e-'^j (184) 

which all yield P[ainin] ~ 0. 

B. The deconflnement phase transition 

Clearly the deconfinement phase transition is related 
to the transition in the effective potential from its IR 
behavior (fT74l) to its UV behavior (|T68| . To illustrate 
this let us approximate the Gribov formula (I151|) by 




(185) 



i.e. we choose the gluon energy as 

w(p) = w„=i (p, A = 0) -I- w„=-i (p, A = 0) . (186) 

This expression is correct in both the IR and UV but cer- 
tainly introduces an error in the mid-momentum regime, 
which influences the deconfinement phase transition. 
With Eq. (dHU) the energy density ([T^ for x{v) = 
becomes 

e"(a,L) = erR(a,i) + e^v(a,i) 



- AJ_ V — 

n=l 



— - {MLf V sin^ (nacrL/2) 

-I cr>0 

(187) 

For SU(2) the only positive root is ct = 1 and this poten- 
tial can be expressed as 

e"(a,i) = \jlfi^)^ fi^) ^ ^\^ - 1)' + <^^ - x) . 

(188) 
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where 



27r2 



(189) 



and the variable x is defined in Eq. (I178p . For small 
temperatures L~^ <C M this potential is negative and 
the system is in the confined phase. As L~^ increases 
the minimum of the potential at x = 2 turns into a 
maximum and the system makes a transition to the de- 
confined phase, see Fig. 2(d) In the relevant interval 



< x < 1 the potential has a single minimum in the con- 
fined phase and two degenerate minima in the deconfined 
phase, which merge to the single minimum at x — ■^ oi 
the confined phase at the phase transition. Starting from 
the deconfined phase the phase transition occurs where 
the three roots of f'{x) — degenerate. This occurs for 
c = h, i.e. for 
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FIG. 5. The energy density (11931) for the gauge group SU(2) 
as a function of the variable x (I189|) for different L~^ in the 
range from 360 to 480 MeV (from bottom to top). The phase 
transition occurs at a temperature of Tc ~ 432 MeV. 



Tr = L- 



■K 



(190) 



With the lattice result M = 880 MeV this corresponds 
to Tc « 485 MeV. 

To exibit the deconfinement phase transition for the 
gauge group SU(3) we cut the potential along a line which 
contains one of the degenerate minima both in the con- 
fined and deconfined phase. A convenient cui|j is y = 0. 
Fig. U shows the y — Q (i.e. ag — 0) cut through the 
potential (I187p . At low temperatures the minimum cor- 
responds to the center symmetric configuration x — -^ 
{y — 0). At a certain critical temperature this minimum 
disappears while a new minimum at a smaller a;-value oc- 
curs, which (for sufficiently high temperature) eventually 
becomes the minimum x = y = for which the Polyakov 
loop equals the trivial center element z = 1. From Fig.|3] 
we find a critical temperature of Tc ~ 597 MeV. 

The critical tempertures obtained above are by far too 
high. Of course, given the approximations used to arrive 
at these values we do not expect a decent description of 
Tc- These approximations are: i) neglect of the curva- 
ture xip) ^-iid ii) approximating Gribov's formula (I185p . 
The use of the correct Gribov formula (J15ip only slightly 
reduces Tc as long as the curvature is neglected as we will 
show now explicitly for the gauge group SU(2). 

Using the Gribov formula (I15ip for a;(fc) but still ne- 
glect the curvature the effective potential (jl62p reads: 



;(a,i):= 



E 



1 



00 



71 — — CK 



^ciPa) 



(191) 



The evaluation of eQ{a,L) has to be done numerically. 
To avoid UV divergences we define 



e^{a,L):^e^{a,L)^e^{a^O,L) 



(192) 



and represent this quantity as 



e^(a, L) = [eg (a, L) - e^y{a, L)] + e^y{a, L) 
= E /<l'p{[^G(pJ-^a=i(p,,A = 0)] 

- [wG(Pa=o) - Wa=-l(p„=o> ^ = 0)] | + 



(a,i) 



(193) 

and calculate the integral numerically while using the 
analytic expression (jl70p in the last term. The resulting 
effective potential is shown in Fig. [5] for the gauge group 
SU(2) for several temperatures L^^ . One finds a critical 
temperature for the deconfinement phase transition of 
Tc « 432 MeV, which is still much too high. 

It is the omission of the curvature, which pushes the 
deconfinement phase transition to higher temperatures. 
This can be seen as follows: In the UV the curvature 
x{p) is suppressed by 1/lnp compared to the gluon en- 
ergy w(p) and is in addition negative. Thus neglecting 
the curvature at most decreases the UV part of the poten- 
tial, which is deconfining. In the deep IR the curvature 
X{p) agrees with w(p) (cf. Eqs. (fT^ and P^ ). Ne- 
glecting here xip) definitely increases the contribution of 
the confining IR potential. So in total, the neglect of the 
curvature x{p) decreases the deconfining part and at the 
same time increases the confining part of the potential. 
Both effects increase the transition temperature. This 
will be confirmed by the numerical results given in the 
next subsection. As we will explicitly see below the crit- 
ical temperatures decreases substantially when the cur- 
vature x{p) is fully included. 



C. The full potential 



Alternative possible cuts are y - 



-i= and y : 



2 



Finally, we calculate the full effective potential (jl62p 
with the curvature x included. We first use the Gribov 
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FIG. 6. The effective potential (fTM)) for gauge group SU(2) 
(a) and the y — {ag — 0) slice for SU(3) (b) as a function 
of X for different L'^ in the range from 260 to 290 MeV in 
figure (a) and from 260 MeV to 300 MeV in figure (b) (from 
bottom to top). We extract the following phase transition 
temperatures: T^ « 267 MeV (SU(2)) and 278 MeV (SU(3)). 




FIG. 7. The effective potential (fT96)) for gauge group SU(2) 
(a) and the y — (ag — 0) slice for SU(3) (b) as a function 
of X for different L~^ in the range from 260 to 290 MeV in 
figure (a) and from 260 MeV to 300 MeV in figure (b) (from 
bottom to top). We extract the following phase transition 
temperatures: n « 269 MeV (SU(2)) and 283 MeV (SU(3)). 



formula (J151I) for uj(p) and in accord with the gap equa- 
tion (jl45p for x{p) the IR expression xmip) (|152p . which 
agrees with the IR behavior of lo {p) . The energy density 
(I162p then becomes 

e(a, L)=J2 [ <J'P [^g(pJ - Xir(pJ] • (194) 

Substracting the value at a = we have 



e(a, L) = e{a, L) — e(a, L ~ 0) 



j(a,L) 



-IR 



ia,L), 



(195) 



where eQ{a,L) and ej^(a, L) are defined in Eqs. (|193p 
and (J174I) . respectively. (The reader is advised to com- 
pare this expression to Eq. (jl87p !) The resulting poten- 
tials is sh own i n Fig. 6(a) for the gauge group SU(2), 
while Fig. |6(b)| shows the as — cut through the SU(3) 
potential for various temperatures L~^. In both cases 
one observes a (deconfinement) phase transition, which, 
however, is first order for SU(2) and second order for 
SU(3). This is because at the deconfinement transition 



the position of the minimum of the potential changes con- 
tinuously for SU(2) but discontinuously for SU(3). From 
these potentials one extracts a critical temperature of 
Tc « 267 MeV and T^ « 277 MeV for SU(2) and SU(3), 
respectively. 

Alternatively we use for xip) the parameterization 
(|154l) fitted to the numerical data and solve with this 
x{p) the gap equation (jl45p for uj{p). To avoid UV di- 
vergencies in the numerical calculations we represent the 
potential e{a,L) (|162p in the form 



eia,L) = ^ d^P w(p„)-Wa=i(p^,A = 0) 

- "(Pct)xir(p<t)J +eQ=i(a,i,A = 0)(a,L) 

-J2 fd'pv{p,)xvy{p.)- (196) 

The first integral can be straightforwardly carried out 
numerically. Using the explicit form of xuv(p) (|153|) and 
v{p) (|155p and the definition (|164p we may write the last 
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FIG. 8. The effective potential (fT94)) for gauge group SU(3) 
as a function of x and y (a) below and (b) above the phase 
transition temperature. In both plots the local minima are 
marked by a contour in the xy-plane. 



FIG. 9. The Polykav loop {-P[a]) evaluated at the minimum 
a = amin of the full effective potential as a function of T/Tc 
for the gauge groups SU(2) (a) and SU(3) (b). 



term as 



Y^ / d-V(P^)xuv(p^) 



E /cr^PxuvlpJ-A-yd^^p^ 



XUv(Pa) 



A)" 



A 



euv(«'^) - ( ]^ j / dtea=i-2t-2n{a,L,X), 

(197) 

where e^y is defined by Eqs. (jB14p of appendix |B] and 
ea{a, L, A) is given by Eq. p66p . The resulting effective 
potentials e(a, L) p96p are shown in Figs. 7(a) and |8] for 
the gauge group SU(2) and SU(3), respectively. Fig. |7(b)| 
shows the y — cut through the SU(3) potential shown 
in Fig. |5] for various temperatures L~^. Figs. |7(a)| and 



7(b) 



and 
for 



. 4^ 

hardly differ from the potentials shown in Fig. 



6(b) 



6(a; 



where the Gribov formula (I15ip was assumed 
(p) and the gap equation solved for x{p) yield- 
ing p52p . The critical temperature found from Fig. [7] 
is Tc « 269 MeV for SU(2) and T^ « 283 MeV for 
SU(3). These critical temperatures only slightly differ 



from the values found in Figs. [S] from Eq. (I194p . Within 
the given accuracy these values agree with the result of 
Ref. [i^l where the variational approach to Yang-Mills 
theory in Coulomb gauge was directly extended to fi- 
nite temperatures, studying the grand canonical ensem- 
ble. The present description has, however, the advantage 
that we do not need an additional variational ansatz for 
the density matrix of the gluons. 

Finally, Fig. [S] (a) and (b) shows the Polyakov loop 
^[o-min] as functions of the temperature obtained from 
the minima Oniin of the full effective potential for the 
gauge group SU(2) and SU(3), respectively. The order 
of the phase transition (second order for SU(2) and first 
order for SU(3)) is manifest in this figures. 



IX. CONCLUSIONS 

In the present paper we have studied the deconfine- 
ment phase transition by investigating Yang-Mills the- 
ory for one compactified spatial dimension in the pres- 
ence of an external background field directed along the 
compactified dimension and living in the Cartan subal- 
gebra. In this formulation the inverse of the length L of 



the compactified dimension represents the temperature. 
The vacuum energy density as a function of the constant 
background field serves as the effective potential of the 
Polyakov loop, the order parameter of confinement. We 
have calculated this potential within a variational ap- 
proach by minimizing the energy density for given (value 
of the) background field and given compactified length 
L, using a Gaussian type ansatz for the vacuum wave 
functional and a background gauge fixing. This resulted 
in a set of coupled DSEs for the gluon and ghost prop- 
agators. We have shown that these equations can be 
reduced to the corresponding DSEs in Coulomb gauge 
at L~^ = and in the absence of the background field. 
However, the momenta along the compactified dimension 
have to be taken at the discrete Matsubara frequencies 
and are shifted by the background field. Using the zero- 
temperature results for the gluon and ghost propagators 
as input, from the effective potential of the Polyakov loop 
we have extracted a critical temperature for the decon- 
finement phase transition of Tc ~ 269 MeV for SU(2) 
and Tc « 283 MeV for SU{3). Within the given accuracy 
this values are consistent with the result of Ref. [l7j . 
where the variational approach to Yang-Mills theory in 
Coulomb gauge was extended to non-zero temperature by 
making a quasi-particle ansatz for the gluon density ma- 
trix and minimizing the free energy with respect to the 
quasi-giuon energy and the Fock basis. The presently 
calculated effective potential has also correctly revealed 
the order of the deconfincment phase transition for SU(2) 
and SU(3). 

In the present approach the deconfinement phase tran- 
sition is entirely determined by the zero-temperature 
propagators, which are defined as vacuum expectation 
values. This shows that the finite-temperature behavior 
of the theory, in particular the dynamics of the decon- 
finement phase transition, must be fully encoded in the 
vacuum wave functional, at least within the present trun- 
cation scheme. 

The results obtained in the present paper are encour- 
aging to extend the present approach to full QCD [15] at 
finite temperature and baryon density. 
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Appendix A: Equations of motion 

In this Appendix we show that the equations of mo- 
tion resulting from the constrained variational approach 
(|i7)l . i.e. the gap equation and the ghost DSE, can be 
consistently solved for propagators of the form assumed 

(ma 



in Eqs. (|112p . (J113p . As in the body of the paper we 
assume a constant background field living in the Cartan 
algebra. We first check that the form (|112p is consis- 
tent with the symmetries of the gluon propagator. The 
symmetry relation (j72p reads in momentum space 



2?.t(p)-^?fc(-p)- 



(Al) 



Using 



tlKp) = tik i-p), (A2) 



Eq. (|74p can be rewritten in momentum space as 

ttt, (p)2?^;?; (p)4?(p) = Kiip) ■ (A3) 

Defining the gluon propagator in the spherical basis by 
Eqs. (jl03p and using (|98)) . (jlOO[) the symmetry relation 
(|A1[) becomes in the spherical basis 



while using (fT05)l Eq. ([M]) reads 



tl,,ip)VtT,ip)tJ,,ip) 



2?fcr(p) 



(A4) 



(A5) 



Using Eq. (jlOSp one easily convinces oneself that the di- 
agonal ansatz (|112p is compatible with Eqs. (jA4p and 

(Ei). 

Next we analyse the equations of motion, i.e. the ghost 
DSE (frfp and the gap equation (|93|) together with the 
expressions ([7^ for the ghost self-energy and ([57)1 for 
the curvature. We start with the ghost DSE ([77]) . 
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the manuscript and useful comments. This work was sup- 
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1. The ghost self-energy 

Rewriting the ghost self-energy ([M]) in momentum 
space 



,bb' t 



-ib'a' 



E(p) = / dq ttt, (-q)C' i^mk^ i<l)Tadk' (P + q)G(p + q)ndi{p) . 



(A6) 



and switching to the spherical basis defined by Eqs. (jl03p 
we find 

j:'^^p)^{-r'+^'fl^,^^f^r'.r 

dqtrA<l)'DZi^Kki^)dUP + ^)G^'^'iP + ^)dJ{p), 

(AT) 



where we have also defined the structure constant T""^ 
^abc jj^ -(-j^g spherical basis 



t,r = e^'''{a\p){b\a){c\T) 

= e'^'^e^^e^ = e^ (e, x e^ 



(AS) 
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Using Eq. (llOip we obtain 

^ par l^ptrr ■ 



(A9) 



Here the totally anti-symmetric tensor epo-r is defined as 
usual, however, for the indices a — — f,0, 1. Since the 
indices p,cr,T,... take the values 0, ±1 the totally anti- 
symmetric tensor Sp^r is non-zero only for 

p+(T + T = 0. (AlO) 

Assuming in Eq. <\A7\i color diagonal propagators (I112[l . 



(|113p we obtain 

^ (P) = ^J- -a'fj.a-l- -a' fiT 

X I dqt';,{q)V^{q)dUp + q)G^'(p + q)d[(p) , (All) 
where we have used 

tfj,(q)if;,.(q)i^,,(q)-tr.(q)- (A12) 

Since (no summation over p, /i). 



T 



-^ —pur ^ar \-^ - 



-pp-cr-^ -PM 



p/j(T| 



(A13) 



the ghost self-energy becomes indeed color diagonal 

E'^^(p)=(5'"^I]'^(p), (A14) 

where 

2 



S'^(P) 



T 



J dqd?ip)tf^{ci)dUp + q)P^(q)G-(p + q) . (A15) 



With the explicit form of d"^ (p) (J107I) we have 



4(P 



<(q) 



'(p)- 



Using furthermore that 

d;:(q)i;:,(q) = o 

we can rewrite the ghost self-energy (|A15[) as 



(A16) 
(A17) 



S'^(P) 



l-L-riial 



X y'dgdr(p)ir,(qK-^(p)P^(q)G-(q + p). (A18) 

The antisymmetric tensor T-Tfia- is non-vanishing only 
for n + a — T = 0. Due to this sum rule the summation 
over T can be explicitly carried out yielding with 



T_ 



1 



(A19) 



- -{ti+a)pa 

the final expression 

S'^(P) = 

- ^ /'d9dr(p)irfe(qK(p)2?''(q)G-+^(p + q) . 

(A20) 



Note this quantity is positive definite since the <i^(p) 
(I107P are purely imaginary. With Eq. (IA14|) the ghost 
DSE (|77| becomes in the spherical basis 



C^' {p)^Gf (p)-S-(p). 



where 



Gf (p) - -d'^(p)d'^(p) = (p - e,aa) 



(A21) 



(A22) 



and S'^(p) is defined by Eq. (TMO)) . This shows that the 
diagonal ansatz, Eqs. (J112p . (|113p for the propagators is 
indeed compatible with the ghost DSE. What remains to 
be done is to show that this ansatz is also compatible with 
the gap equation (|93l) . For this purpose we investigate 
first the curvature x(q): which enters the gap equation. 



The curvature 



Fourier transforming the curvature (j87p we find 



xtiiP) = ltii,ip)ti^:{^p) 



J dqtr (fa'dk' (p + q)G(p + q)fvdr (q)G(q) 



(A23) 



Rewriting this expression in the spherical basis (jl03p and 
assuming a color diagonal ghost propagator (I113P we ob- 
tain 

Xkl (P) ^ ^-'- -rpa--'- -rva 

X / dgi;:,,(p)4,(q + p)4(q)ir.;(p)G'^(q)G^(q + p). 

(A24) 



With (IA13P we find that the curvature is also color diag- 
onal in the spherical basis 



xtiiP)-S^''xUp)^ 



(A25) 



Using (|A16p and the sum rule for the indices of TLr^o- the 
summation over r can be trivially carried out, yielding 

xUp) = 

-J dp t^,,, {p)dl, {q)df, (q)ir,, (p)G- (q)G^+- (q + p) . 

(A26) 

Consider now the Lorentz structure of XkiiP)- ^Y its 
definition as derivative with respect to the (transversal) 
gauge field this quantity has to be transverse as well. 
Indeed from the explicit expression (|A26p it is seen that 
it satisfies the transversality condition 



c.(p)x;:,(p)irn(p)=x^„(p) 



(A27) 
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From this we can conclude that it has the fornri 

X^,(p) = ^t^,(p)x"(p), (A28) 

where the scalar curvature x^(p) can be obtained as 



^"(p) = d3i*r.(p)x;:,(p). 



(A29) 



Here d is the number of spatial dimensions and we have 
used 



£i;,(p) = d-i. 



(A30) 



Inserting here the explicit expression (IA26|) and using 
*n(p) — *fz'(p) ^c obtain 



X''(P) = - 



2(rf-l) 

J2 I dqdimKip)dri<DG^Ci)G'^+^p + q) . 

(A31) 



Eqs. (JA25I) and (|A28p show that in the spherical basis 
the curvature x{q) is not only color diagonal but has 
also the same Lorentz structure as assumed for the gluon 
propagator pi2p . 



structure of the gluon propagator assumed in (|112p . 



which resulted in the Lorenz structure (jA28p of the cur- 
vature. Thus the ansatze (J112p . (J113p for the propagators 
are also compatible with the gap equation. We have thus 
shown that the equation of motions resulting from the 
constraint variational principle (H)^ — )• min, (A)^ = a 
in the background gauge (j67p allow for color diagonal so- 
lutions. For such solutions we find for the energy density 
([M)) at the stationary point 



(HK + Hi)^ 



5VE 



y'dgQp'^(q)-i-X^(q)) , (A34) 



where V is the spatial volume and 'D°'{p) and x'^ip) are 
still Lorentz matrices. Carrying out the trace over the 
Lorentz indices thereby using Eqs. (J112p and (|A28p we 
obtain 

(Hk + H^)^ 

= (d-l)V^^y'dg(^^P'^(q)-i-x'^(q)) , (A35) 

where 'D'^ip) and x'^(p) are now Lorentz scalars. (Note 
by Eq. (jA28p the Lorentz tensor xliiP) ^^d the Lorentz 
scalar x'^(p) differ by a factor of g^.) 



3. The gap equation 

Consider now the gap equation (j93p . which after 
Fourier transformation becomes 



V-\p)V-\p) 



£(p)d(p)d(p)£(p) + g\ip)x{p) ■ (A32) 



9 



This is still a matrix equation both in the adjoint color 
and Lorentz indices. Rewriting this equation in the 
spherical basis defined by Ea. (jl03p and using Eqs. (|104p 
and (IA25P we observe that the gluon propagator I?(p) 
has indeed to be color diagonal and obtain 



= (-r(p)d-(p)d-(p)t-(p)) + .gV(p)x^(p) . (A33) 

Here 2?'^(p)^^ and x'^(p)i like ^"^(p), are still Lorentz 
matrices. This equation is compatible with the Lorentz 



Appendix B: The energy density 

Below we calculate the contribution of the generic 
gluon energy (|164p to the energy density (|165l) . Rep- 
resenting 

( VpFTa) " = (pi + {p„ - craf + a) "^' (Bl) 
by a proper-time integral^ n 

ipl + ^r'' = 777^ drr-i-te-(^^+^) (B2) 

i [-2) "'O 

we find for the energy density (|165p 



ea(a,i, A) 



4- r(-f) 



J2j:J2 d^/ drr-i-te-^[-+('="-'^'^)']e-^^ 



a n 



(B3) 



Here we have carried out the integral over the azimutal 
angle of p^ and set x — p\. The integral over x can 



The transversal projector i^;(p) will explicitly show up after the 

angular integral in llA26l l is carried out. Furthermore, the factor " We assume here that a is such that the proper-time integral 

g2 has been explicitly included so that x^(g ) redu ces for ^ = exists, otherwise f^ dr has to be replaced by f^^2 dr, where A 

to the curvature defined in Ref. [^, see Sect. IVICI is a UV cut-off. 
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be carried out yielding a factor r ^. By raeans of the 
Poisson resummation formula one derives the relation 



1 1 r°° 






2mr 

17' 

(B4) 

Using this relation in Eq. (|B3I) and performing the in- 
tegration over z we obtain 



60(0, L, A) 



A-f 



/fl — a 



(47r)3/2r(-f) 



cr n— — 00 



Obviously the proper-time integral is not well-defined for 
A = and n = 0. Fortunately, this term drops out by 
subtracting the energy density at vanishing background 

I 



field 



ea{a, L, A) := ea{a, L, A) — ea{a = 0, L, A) 
Afi-" 2 



(4^)3/2 r(-f) 

X V" V" [cos(no-aL) - 1] / drr^^"^c^^(^) ~^'^. 
<T „=i ■^0 

(B6) 

With a change of variable s = {nL/2) /r in the proper- 
time integral 



r dTT~i'^e-^^^)\-^^ . (B5) / 

Jo ■'^ 



drr 



5 a lfnL\2 -., 



2 

7l 



3+a 



/ dss2+-e--i — j (B7) 



we find for the energy density (jB6p after some trivial 
algebraic manipulations 



e„(a,i,A) = -4(4^)-3/2 



A/i 



r(-f) 



2^ 2^ \ „T. 



a+3 



sm^ f n I / dss2 + Se '^e 



-)Jo 



)V^r 



(B8) 



For A = the remaining integral can be expressed by 
Euler's Gamma-function 



ea{a,L, A = 0) 



EvL 



4Af 1-" r (I -H f ) 

"(47r)3/2 r(-f) 

E;3Tirsin'(n<TaL/2). (B9) 



3+a cxD 



n=l 



Obviously, the vanishing roots a — (0, 0) do not con- 
tribute to the potential difference ea{a) — Caia — 0). 
Furthermore, the negative roots give the same contribu- 
tion as the positive ones. We can therefore restrict the 
summation over a to the positive roots ct > and take 
care of the contribution of the negative roots by an ad- 
ditional factor 2. This yields 



60(0, i, A = 0) = —8 



Afi-" r(|-Hf) /2 



3+a 



(47r)3/2 r(-f) \L^ 
EE;3T^«"^'("-^V2). (BIO) 



a->On=l 



For X ^ the remaining proper-time integral in Eq. 
cannot be taken analytically but can be expressed by a 



r 



modified Bessel function K^{z) 



dss ^ *e ^ exp 



m) 



yielding 
ea(a,L,A) 

00 

xEE 



= 2(^y=] KtinLVX), (Bll) 



Ml 



'2n/a' 
nL 



(47r)3/2r(-f) 

2^2 



sin^ (crai/2) K_^_3{nLV\) 



(7 n=l 

(B12) 

As in the previous case A = 0, the summation over a can 
be restricted to cr > by including a factor of 2. 

To eliminate the UV divergencies from the energy den- 
sity e(a, L) (|162l) we also need to isolate the contribution 
from the UV behavior of x- The UV form of x is given 
by Eq. (|T53| . Using 

f2 \* 



1 






dt 



Af2 



p2+A 



, p2 ^ A > M2 (B13) 



and the definition (|164p of oJq (p, Af , A) we have 

/•oo 

Xuv(p) = / dia;a=i_2t(p,A). (B14) 

Jo 
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With this relation the (negative of the) energy contribu- 
tion generated by the UV-behavior of x reads 

el^{a,L):^Y.j ^P^^^^P-) (^15) 



where ea{a, L, X) is defined in Eq. (jB3p . Subtracting the a = contribution and using Eq. (|166p we obtain 

M2* /Va\'"*^^/2\2-* . ,/ L 



-uv 



(a,L)-e^v(a = 0'^) = -- 



df- 



(4^)3/2 7o r{t-l)\L 
and after restricting the summation over a to positive roots 



'UV 



{a,L)-e^y{a^O,L) ^ 



EE 




2-t 



EE 

cr>On=l 



oo / n\ 2 — t 



sin^ I ncra— j Xf_2 (nLvAj 

(B16) 



^2 ' ncra| ) i^t_2 (nlVX 



(B17) 
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